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We investigate the influence that adding a new coupling has on the linear stability of the synchronous state in 
coupled oscillators networks. Using a simple model we show that, depending on its location, the new coupling 
can lead to enhanced or reduced stability. We extend these results to electric power grids where a new line 
can lead to four different scenarios corresponding to enhanced or reduced grid stability as well as increased or 
decreased power flows. Our analysis shows that the Braess paradox may occur in any complex coupled system, 
where the synchronous state may be weakened and sometimes even destroyed by additional couplings. 
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I. INTRODUCTION. 

Collective synchrony is an omnipresent phenomenon in 
systems of coupled oscillators [1,2]. It arises when the cou¬ 
pling between individual oscillators becomes strong enough 
that it overcomes the tendency of oscillators to swing at their 
natural frequencies. Simplified models such as the Kuramoto 
model [1,3] allow to describe a plethora of non linear phe¬ 
nomena involving collective synchrony in Josephson junction 
arrays [4], biological systems [5, 6], crowd dynamics [7], cou¬ 
pled neural networks [8], chemical reactions [9] and elec¬ 
tric power grids [10-13] to name but a few. Quite naturally, 
one expects that adding couplings between initially uncou¬ 
pled pairs of oscillators generically favors synchrony. This 
is however not always the case, Nishikawa and Motter pro¬ 
vided analytical conditions for systems of coupled oscillators 
to be synchronizable over a larger parameter range [14]. Ref¬ 
erences [15, 16] found numerically that adding a new cou¬ 
pling in an initially synchronous system sometimes destroys 
synchrony. This unexpected scenario is the electrical analog 
of the Braess paradox, first discussed in the context of traffic 
networks [17, 18], where building new roads sometimes in¬ 
creases traffic congestions. Similar counterintuitive observa¬ 
tions were reported for simple mechanical systems and uncon¬ 
trolled electric circuits [19, 20]. One purpose of the present 
manuscript is to present a more systematic analytical treat¬ 
ment of the Braess paradox in coupled oscillator systems. 
While our focus is on electric power grids, our theory also 
applies to other oscillator networks described by similar mod¬ 
els. 

The operational state of AC power grids requires syn¬ 
chrony of thousands of rotating machines of widely varying 
sizes, millions of electric and electronic devices and compo¬ 
nents, over several voltage levels intercoupled by frequency¬ 
preserving transformers [21]. Nowadays, power grids are 
maintained in a synchronous state at their rated frequency (50 
or 60 Hz) by active control of power generators. An inbal¬ 
ance between production and consumption results in a vari¬ 
ation of the operating frequency but not necessarily in the 
loss of synchrony. The latter may arise if frequency variations 
exceed safety margins that require to disconnect parts of the 
network. The standard operational protocol is crucially chal¬ 
lenged by the current rise of weakly controllable renewable 
energy sources. Maintaining the operational state and guaran¬ 


teeing the safe distribution of power under these changing cir¬ 
cumstances requires power grid upgrades, in particular the ad¬ 
dition of new transmission lines. This is however both costly 
and not always well accepted socially. It is therefore crucial 
to upgrade grids efficiently, adding as few lines as possible to 
ensure a more stable and safer grid operation. Understanding 
the Braess paradox in electric power systems is therefore key 
to optimize the grid of tomorrow. 

Improvement in grid operation after a line addition can be 
quantified for instance by (i) the new power flows on initially 
strongly loaded lines, this measure of power re-routing is re¬ 
lated to line outage distribution factors used in electrical engi¬ 
neering [22], (ii) the linear stability as measured by the Lya¬ 
punov spectrum [12, 23-26] of the upgraded grid and (iii) the 
size of the basin of attraction of the synchronous state in the 
associated parameter space [27]. Further criteria include N-1 
feasibility and voltage stability [21]. In this work we investi¬ 
gate the impact of line addition on grid operation along points 
(i) and (ii) in a purely reactive power grid. We illustrate ana¬ 
lytically on a simple chain network how the perturbative addi¬ 
tion of a line, which modifies the grid topology by creating a 
loop, affects the power load of the electrical connections and 
the linear stability. We classify the impact of the new line into 
one of four different scenarios - depending on whether lin¬ 
ear stability is improved or not, and whether strongly loaded 
power lines are relieved or not. Out of these four possible 
scenarios, three are different manifestations of the electrical 
Braess paradox, where (I) already strongly loaded lines be¬ 
come even more strongly loaded, (II) network stability is re¬ 
duced or (III) both. We furthermore show how these three sce¬ 
narios for the Braess paradox also occur in a complex network 
having the topology of the British electric transmission grid. 
Our analytical calculation contributes to the understanding of 
the Braess paradox in electrical systems. We conjecture that 
the paradox is generic and may occur in any system of coupled 
oscillators with reduced connectivity. 


II. THE CHAIN MODEL. 

We consider an AC electric power system in the form of a 
chain connecting N + I nodes [see Fig. 1]. A unique gener¬ 
ator (labeled i = 0) is located at one end of the branch while 
the remaining N nodes (labeled i - 1,..., A) are all loads. A 
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necessary condition for the system to be in steady state is that 
the total injected power at the generator is equal to the total 
power consumed by the loads. For a power injection Pq > 0 
at the generator and F, < 0 at the loads, in arbitrary units, 
this amounts to 
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FIG. 1. The chain model. A single generator (square) injects a power 
Po > 0 which is consumed by N loads (circles) each consuming a 
power of Pi < 0 in arbitrary units. The lines have capacity > 
f’o“2/=i IFj, except the newly added line (dotted) which has capacity 

S. 

age transmission grids, the line admittance is dominated by 
the susceptance (its imaginary part). Accordingly we neglect 
ohmic effects and assume that each node has a constant inter¬ 
nal voltage magnitude |y,j. Under these approximations the 
active power flow equations read [15, 21, 28-31] 
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III. IMPACT OF LINE ADDITION ON POWER FLOWS. 


To leading order in 6/K, the power flowing through the ad¬ 
ditional line is o ^ sin q [32]. The sign of and thus 
the direction of the power flow, changes as a function of d. 
The power flowing through the (0,1) line between the gener¬ 
ator and the first node goes from Pyo = Ao,i sin0i_o = -Po 
to Pi,o = ^ 1,0 sin §10 ~ -Po - <5 sin 0^/ 0 once the new line is 
added. As long as sin Od o < 0, the new line lowers the load 
on all the lines {i, / + 1) for i = 0,... ,d - 1. However, when 
sin 0d,o > 0, we face the counterintuitive situation where the 
new line transmits power back from node d to the generator, 
thereby increasing the load on all the lines between the gen¬ 
erator and node d in the original network. This is an electric 
manifestation of the Braess paradox [17, 18] and its occur¬ 
rence is due to the nonlinear nature of the power flow Eqs. (1). 
In the case of our simple model, which of these two scenarios 
takes place depends only on the value of Odfi- 


IV. LINEAR STABILITY. 


,ysin(0j-0i) ie{0,l,...,N} , (1) 

H 

where j ~ i indicates that the sum over j spans the nodes 
connected to the node, Kjj - Kji = Bj j\Vi\\Vj\ denotes the 
maximum power capacity of line {i, j) having susceptance Bjj 
and 0, is the voltage angle (with respect to the curi'ent) at node 
i. Since P* = Pq - \Pi\ units of power are transmitted 
from node i to (-H1, the line capacities must satisfy A,+i, > P* 
for Eqs. (1) to have a solution. Solving Eqs. (1) for the angles 
yields 

0i+i,i = Oi+i - 0/ = - arcsin [p*/A;+i,/] i e (0,1,..., A - 1], 

( 2 ) 

so that 0,+y,' 6 [-71/2,0]. 

We next add to this topology a line of capacity 6 be¬ 
tween the generator and the t/* load, d e (1,..., A]. When 
6 <si Ki+id for all /, the perturbed solution [0,] remains close to 
the unperturbed one, i.e. 0,+i_,' « 0,+i_; -H with |e,+i,,| « 1. 
Solving for the e,+i_,’s, one obtains explicitly the order cor¬ 
rection to the unperturbed power flow solution as 

e<+l,; = j COS0;+1,; (3) 

[0 i>d. 

Clearly, e,+i, = 0 for / > li since the power flowing through 
the lines connecting nodes i and i + 1 for i > d is left un¬ 
changed. Since the above result is valid to order in d/K, 
all angle differences entering Eq. (3) are differences of the un¬ 
perturbed angles [0,]. In particular, the difference between the 
voltage angles of the nodes which are connected by the new 
line is 0d,o = Zto “ Zfjo ^^sin [p*/A,+i_,] . Be¬ 

low we investigate how the power flowing through the lines 
changes when adding the new line. 


The solutions of the power flow Eqs. (1) describe the oper¬ 
ating stationary state of the power grid at a given time. Upon 
changing conditions, such as variations of the power injected 
and consumed, the angles’ dynamics in this transient stability 
problem is governed by the swing equations [13, 21] 

m + Dfii ^Pi + Y^ Kij sin (0,. - 0,) , (4) 

i~i 

which describe the power balance at nodes with rotating ma¬ 
chines as generators or loads. Without inertia /, s 0, Eqs. (4) 
reduce to a Kuramoto-like model [1, 3, 33], with reduced 
connectivity. Linear stability in the Kuramoto and similar 
models with reduced connectivity has been investigated in 
Refs. [30, 34-36], which derived bounds on the exponential 
rate of return to the stationary state. 

For /,■ + 0, linearizing Eq. (4) around a stationary solution 
©(f) - 6 + 59{t) yields the eigenvalue equation 

MSe = A [Adiag(I) + diag(D)] 86 , (5) 

where A's e C are the Lyapunov exponents of the dynamics 
governed by Eq. (4), diag(I) and diag(D) are diagonal ma¬ 
trices, diag(I),',' = /, and diag(D)„ = D,. Einally, M is the 
stability matrix defined by Mjj - K,j cos 0jj if i and j are 
connected, M„ = - and zero otherwise [12, 23]. The 

stationary solution is linearly stable if the largest nonzero Lya¬ 
punov exponent is negative and unstable otherwise. We next 
show that the system is stable if M is negative semidefinite 
and that the loss of stability occurs when the largest nonzero 
eigenvalue of M becomes positive. This justifies our use of 
the spectrum of M as measure of stability, keeping in mind 
that time scales, e.g. for restoring synchrony may depend on 
/, and Di. 

In the case of homogeneous inertia /, s I and damping co¬ 
efficients Di = D, Eq. (5) is diagonalized by the eigenvectors 
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of M and the Lyapunov exponents are simply given by 


V. LINEAR STABILITY FOR THE CHAIN MODEL 


A 


± 
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( 6 ) 


In the case of the chain model prior to line addition, M is a 
(A + 1) X (A + 1) tridiagonal, symmetric matrix. 


where p - Dj I and is one of the eigenvalues of the stability 
matrix M. In the inhomogeneous case, projecting Eq. (5) onto 
SO gives. 


O^A^a + Ab-c, (7) 

where we introduced the shorthand notation for the overlaps 
a - <50diag(I)<50, b = ^0diag(D)^0 and c = S6MS9, with 
coefficients a,b >0 since inertia and damping coefficients are 
positive quantities (i.e. Z),-,/, > 0 Vi). The Lyapunov expo¬ 
nents then take the form 

A* = — (^-b + + 4acj . (8) 

In both Eqs. (6) and (8), Re[A ] is always negative and stabil¬ 
ity depends on the sign of Re[A^]. 

In the homogenous case, Eq. (6) makes it clear that lin¬ 
ear stability is determined uniquely by the spectrum of M: 
Re[Ag] and Aa become positive simultaneously. Since M is 
real and symmetric, all da’s are real. Thus, the necessary 
condition for the system to be stable is that all Aa are nega¬ 
tive. Furthermore, Re[A^] is negative and decreasing as Aa 
decreases in the interval 0 > do > -/y6^/4, while Re[Ag] sat¬ 
urates at -jS/2 when do is decreased further. Thus, not too far 
from loss of stability, the spectrum of M is as good a measure 
of increase/decrease of stability as the true Lyapunov spec¬ 
trum. 

We extend the approach of Ref [26], which deals with inho¬ 
mogeneous damping but identical inertia, to the case of inho- 
mogenous inertia and damping. When M is negative semidef- 
inite, the coefficient c in Eq. (8) is negative. Thus, given that 
a > 0, Re[A^] is negative and the solution is linearly stable. 
Furthermore, the cancellation of the Lyapunov exponent oc¬ 
curs only when c vanishes. For this to take place, SO must be 
proportional to one of the eigenvectors of M associated to a 
zero eigenvalue. This shows that a stationary solution of the 
dynamic system (4) is linearly stable as long as M is negative 
semidefinite and the loss of stability occurs when the largest 
nonzero eigenvalue of M vanishes. 

We therefore take from now on the spectrum of M as a mea¬ 
sure for increased (Aa decreases) or decreased stability (Aa 
increases). Because inertia does not influence the stationary 
state, taking this latter criterion allows to make more gen¬ 
eral statements regarding stability, however one needs to keep 
in mind that /, and D, may in principle affect stability in a 
nontrivial way. We defer investigations of this issue to future 
work. 

Having discussed the role of the spectrum of the stability 
matrix on the Lyapunov exponents we next investigate how 
stability is affected as a line is added to an initially stable net¬ 
work. 


0 

0 : 

0 

-Cjvjv~i 

■CjV,N~l Cv,iV-l , 

(9) 

_j e [-7r/2,0], M is 
diagonally dominant [37] with only negative diagonal ele¬ 
ments and positive subdiagonal elements. It thus belongs 
to the family of Jacobi matrices [38], in particular M has 
distinct eigenvalues. By Gershgorin circle theorem [39], it 
is negative semi-definite. Furthermore, = (!,...,!) is 
the eigenvector associated to the eigenvalue which vanishes 
by rotational invariance. We order the eigenvalues of M as 
A] = 0 > 42 > ... > Ta^+i. The semi-negativity of the stability 
matrix indicates that the power flow solution of the original 
network topology is stable against small perturbations. Since 
the largest eigenvalue Ai vanishes, stability is determined by 
42. Next, we therefore calculate the leading order correction 
to 42 resulting from the line addition. 

The stability matrix M after the new line has been added, 
has a very similar structure to M except that, first, the angles 
entering in M are the 0,’s and second, the new line modifies 
the following matrix elements: j = -Ci_o - Cd,(h Md+\,\ - 

Mi^d+\ - Cd,o and Md+i,d+i - -Cd-\,d - Cd+i,d - Cd,Q, where 
Cjj = Kij cos Ojj and Kdfi - Ko^d - S. Using Eq. (3), we 
express M as M - M + AM + 0[(5IK)^], where AM is the 
leading order correction to the stability matrix. 


M = - 


Ci.o -Ci,o 0 

-Cifl Ci_o + C2,l -C2,l 

0 

: 0 

V 0 ... 0 

■ ■ = fT ■ ■ rrm 


AM - 6 sin Odfi 


AIM(^+l)x(rf+i) 0(d+l)x(N-d) 
0(N-d)x{d-hl) 0(N-d)x(N~d) 


with AIM defined as 


( 10 ) 


'-CTd,o-Ti,o Ti,o 0 .. 

Ti.o -Tifi-T2,i T2 ,i 0 

0 


CTdfi 

0 


i 0 ■■■ ■■■ Td,d-1 

, CTdfi ... 0 Td4-i -CTdfl - Td4-\ , 

( 11 ) 

where we introduced the notations Tj+ij = tanfl,+i/ and 

CTdfi = cotOdfi. 

Let e be defined by = 42U^^*. Then, the 

leading order correction to 42 is given by A42 = AMvP'\ 
If the sign of A42 is negative (positive), then, to U‘ order in 6, 
the stability of the power flow solution is enhanced (reduced). 
Below we discuss how sgn(A42) changes as a function of the 
position d of the additional connection. To achieve this, we 
distinguish the two cases tan Odfi ^ 0 and tan Odfi ^ 0- 

When tan Odfi ^ 0, the matrix AIM is diagonal dominant [37] 
(since 6 [-7r/2, 0] we have tan0/+i,,' < 0) with a strictly 
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positive diagonal. Hence AM (10) is either semi-positive or 
semi-negative definite depending exclusively on the sign of 
sinflj/ o- In both cases the sign of AA 2 = is well 

defined regardless of and we have 


( A/I 2 > 0 for 6^0 e [njl, tt] ^ reduced stability, 

A/I 2 < 0 for 6 do ^ [-^/2,0] enhanced stability. 

( 12 ) 

When tan 0^0 ^ 0, AM is no longer diagonal dominant 
and it is not possible to determine the sign of AA 2 as directly 
as before. Instead we use = (M[f \ ■ • ■, mJv') to compute 
u(2)'"aMu^^* explicitly. 


AA 2 


-6 cos ddfi 

d-l 

+ tan 0 dfi Yj (“P’ “ 
1=0 


( 13 ) 


The one dimensional nature of the model allows to express the 
difference as the telescopic sum “ 

“(f* - u^f'- Using this identity we rewrite the term (Mq^^ - 
entering Eq. (13) as 


d-l 


d-l 


/=0 J /=0 


d-l 


.2 

/=0 /■>/ 

(14) 

Substituting Eq. (14) in Eq. (13) finally yields 


A/I2 = -6cos0d,o 


d-l 


2 Z - “f+\) 


/=0 j>i 


d-l 


- tan Odfi - M-+i)^ (tan 6>,+i,; -H cot Odfi) 


i '=0 


( 15 ) 

Because M is a Jacobi matrix, it can be shown (See Ap¬ 
pendix A) that the components of its eigenvector are 
monotonously ordered. Thus (m® - )(m® - > 0 

and the first term in Eq. (15) is positive. Eurthermore, for 
tan 0^0 S 0, it is possible to establish a sufficient condi¬ 
tion on 0 dfl according to which the sign of AA 2 is known. 
Since all Oi+ij belong to [-;7r/2,0] and given the monotonic¬ 
ity of the tangent function over this interval, if (tan0™“ -H 
cotfljo) S 0 where = min,£|o_..._d-i)®i+i,i then we also have 
(tan0,+i ,■ + cot0dfi) > 0 for / € {0,1,... ,(i - 1). When this 
is the case we conclude that AA 2 oc - cos Odfl. Thus, when 
tan 0 dfi > 0, the sign of AA 2 is directly given by -sgn(cos 0 dfi) 
if tanflj o ^ “ cot0™’". Inspecting Eq. (2) one sees that 0™” is 
realized on the most loaded line prior to the network upgrade, 
that is on the line having the largest ratio of transmitted power 
over available capacity. Let {q+l,q} denote the line minimiz¬ 
ing min,£|o_..._<^_i)0/+i,; (i.e. the line where P*fKj+ij is maxi¬ 
mal), then the condition tan 0 dfl < - cot 0™’" can be rewritten 
as 


i^n0d,o< ( 16 ) 


This defines a critical angle a = arctan -y ^^+1 ^ - Pq^/Pq \ ^ 
[0,7r/2], such that 


( AA 2 < 0 for 0 dfl 6 [0, a] ^ enhanced stability, 

A/I 2 >0 for 0£/o 6 - 7 T + a] ^ reduced stability. 

(17) 

The size of the region [a, 7r/2] 1J[“^ + “^/2], where the 

evolution of the stability remains undetermined, vanishes as 
P*IKg+i q when P*IKq+i ^ > 0 since in this limit a 7 t /2 - 

Pq/Kq^i,q- 

These results are summarized in Eig. 2. When the an¬ 
gle difference between the newly connected nodes satisfies 
Sd.o 6 [-7r/2,0], the additional line reduces the load on the 
most loaded line in the loop (i.e. line {q +l,q}) and the stabil¬ 
ity of the power flow solution is enhanced. This is what one 
generally expects of line addition. Line addition can however 
worsen the operating conditions of the network and our the¬ 
ory highlights three different Braess scenarios how this may 
happen - by either increasing the load, by reducing I/I 2 I, or 
both. The worst case scenario occurs when 0dfi e [7r/2, tt]. 
Then, the additional line increases the power load on the most 
loaded line and the stability of the new solution is decreased. 
Paradoxical situations occur when Odfi e [0, a] (respectively 
Sd,Q G [-tt, -tt + a]) as the load on the most loaded line in¬ 
creases (decreases) while the linear stability is enhanced (de¬ 
creased). These three outcomes are three different manifesta¬ 
tions of Braess’s paradox in electric power transmission. We 
note that the chain model results also apply to the case of line 
addition in radial (tree like) networks as long as the new line 
connects two nodes on the same branch. 



FIG. 2. (Color online) Impact of perturbative line addition on the 
linear stability of the power flow solution (green region, enhanced 
stability; yellow region, reduced stability) and on the load of the 
transmission lines (top quadrants, increased load; bottom quadrants, 
decreased loads) as a function of the value of 0j,o- 


VI. EXTENSION TO COMPLEX NETWORKS. 

To show how the mechanisms described above can lead 
to the loss of synchrony in more complex networks, we 
consider the electric power transmission grid discussed in 
Refs. [15, 16]. It has the same topology as the UK transmis¬ 
sion network, and we take the same distribution of loads, gen¬ 
erators and line capacities as in Ref. [15]. The general struc- 
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ture of the grid is that of a northern, importing zone connected 
to a southern, exporting zone via only two lines which are al¬ 
most at full capacity [see Fig. 3 (left panel)]. It is obviously 
desirable to relieve these lines by adding another south-north 
transmission line. 



FIG. 3. (Color online) Left: UK transmission grid, 10 generators 
with power P = 11 (squares), 110 loads with P = —\ (circles) and 
uniform line capacity K = 13. Power flows are represented by arrows 
and their magnitude is color coded. The dashed lines (a) - (c) rep¬ 
resent three different line additions considered and the solid line de¬ 
notes the network partition into northern and southern zones. Right: 
Plot of the difference in power flows between the solutions after and 
before the addition of line (a) of capacity S = 1.5. Arrow heads are 
drawn only for power flow differences larger than 0.01. 

The situation is in a way similar to our simple model, where 
the south plays the role of the generator and the north that of 
the loads. This is however only an analogy since the elon¬ 
gated UK grid is a meshed network and not a ID model as 
considered above. In the case of a generic network, Eq. (15) 
becomes 


AA 2 = d sin ffafi 


fj [uf - uff - (m® - cot 

{i.j) 


(18) 

with fij = K,jsinfl,j 2 z> 2 (“® - “‘I 

where a and fi are the nodes connected by the new line, (/, j) 
indicates the sum over all pairs of connected neighbors in the 
original network and is the /* eigenvector of the stabil¬ 
ity matrix (see Appendix B). After the upgrade, the power 
flowing through the line connecting nodes i and j becomes 
Pi j - Kj jSinOj j + 6 sin Oafi fijCotOjj. Following our work, 
a similar expression for the change in power flows resulting 
from the variation of the capacity of one line was used in 
Ref [40] to investigate the effect of line failures. 

In what follows we present examples of additions of new 
lines between the north and the south. Each illustrates the 
realization of one of the electric Braess paradoxes discussed 
above. We first add the dashed line {a) [Eig. 3], between two 
nodes having the angle difference flNorth “ ^South = Sn.S ~ 
0 . 97 T 6 [7r/2,7r]. For this choice. Fig. 2 predicts counter¬ 


intuitively, that power will flow from the north to the south 
through the new connection. This is numerically verified in 
Fig. 3 (right panel), which shows that adding the new line in- 



FIG. 4. (Color online) Lyapunov exponent (dashed) and power flow¬ 
ing through the lines connecting the north and south areas as a func¬ 
tion of the capacity of the additional line d. Each of the panels refers 
to one of the line additions represented in Fig. 3 (labels correspond 
to the labels of the additional lines in Fig. 3) and illustrates one of the 
three Braess scenarios identified in this work. Interestingly, panel c) 
shows the coexistence of two different stable solutions for d > 5.98. 

creases even further the load on the two original lines connect¬ 
ing the two zones - the power flow in the new line goes in the 
wrong direction. The effect is quantified in Fig. 4 a) as a func¬ 
tion of the capacity of the new line. Both the loads on the orig¬ 
inal connection lines and the Lyapunov exponent A 2 increase 
as a function of 6 . Going beyond the validity of our perturba¬ 
tive approach, synchrony is lost in the interval h 6 [2.1,5.2] 
(the gray region in Fig. 4 a), and is recovered at larger values 
of (5, where the stable operating state strongly differs from the 
unperturbed one. Synchrony is lost for 6 y 2.1 and 6 \ 5.2 
as the power flow solutions become unstable (A 2 0), simi¬ 

larly to results reported in Ref. [26]. 

The added line labeled by {b) [see Fig. 3] is chosen to con¬ 
nect two nodes such that (^N,S ~ -0 .971 e [-TT,-7r/2]. As can 
be seen in Fig. 4 b) power is flowing from the south to the 
north along this new line. Despite the associated reduction of 
the power flow on the two original lines, the Lyapunov expo¬ 
nent increases as predicted in Fig. 2. For larger added capac¬ 
ity, however, A 2 reaches a maximum, then starts to decrease, 
and synchrony is never lost. This observation can be under¬ 
stood qualitatively in terms of our simple model: as the capac¬ 
ity of the new connection increases, the difference > orig¬ 
inally in the 3“* quadrant increases until eventually it reaches 
the 4* trigonometric quadrant for which the correction to A 2 
is expected to become negative. 

We Anally add line (c) [see Fig. 3] between two nodes with 
Sn.s ~ 0 . 377 . This time power flows through the new line 
from the north to the south. The loads on the original con¬ 
nections between the two zones, which were already close 
to saturation, increase further. Quite interestingly, the linear 
stability of the solution is improved, AA 2 < 0, despite this 
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load increase. The solution followed as 6 is raised from 0 re¬ 
mains linearly stable in the whole capacity range investigated 
6 6 [0,13]. When the capacity of the additional line reaches 
5.98 units of power, the numerical simulations also converge 
to another stable solution of the power flow equations [see 
Fig. 4 c)]. The behavior of A 2 indicates that the new, large-^ 
solution becomes unstable (A 2 = 0) when 6 \ 5.98. The 
regime h > 5.98 is an example of coexistence of multiple sta¬ 
ble power flow solutions [31, 41^4]. 
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FIG. 5. (Color online) Lyapunov exponent (dashed) and power flow¬ 
ing though the lines connecting the north and south areas as a func¬ 
tion of the capacity of the additional line 6, for the UK transmission 
grid in the case of line capacities uniformly distributed in the in¬ 
terval [9.75,16.25]. Each panel refers to one of the line additions 
represented in Fig. 3 and illustrates one of the three Braess scenarios 
identified in this work. 

To assess the robustness of the perturbation theory results 
we repeat the numerical simulations for the UK transmission 
grid including line capacity variations of +25% with respect 
to the uniform case, K = 13, presented above. We take line 
capacities uniformly distributed in the interval [9.75,16.25] 
keeping the most loaded line crossing the north-south border 
atK- 13, and consider the same network upgrades discussed 
earlier. The numerical results presented in Fig. 5 are very sim¬ 
ilar to those of Fig. 4, indicating that the three different Braess 
scenarios identified for the uniform line capacity case are ro¬ 
bust with respect to the significant capacity variations consid¬ 
ered. 


VII. CONCLUSION. 

We classified the impact of a line addition in an AC power 
grid into four possible scenarios depending on the change in 
linear stability of the synchronous solution and on the change 
in power load on the lines. For the chain model, we showed 
that the effect of such network upgrades depends uniquely on 
the voltage angle difference between the nodes connected by 
the new line. This classification is summarized in Fig. 2, and 


we showed that it can be extended to meshed networks hav¬ 
ing the topology of transmission grids. In this case it is how¬ 
ever less straightforward to predict from the unperturbed op¬ 
erational state which of the four scenarios will be realized. 

We think that our theory has significantly deepened our un¬ 
derstanding of the Braess paradox in electric power systems. 
More generally, it is based on rather generic models, which 
suggests that Braess paradoxes, in the form of weaker stabil¬ 
ity of the synchronous state after coupling addition, are ubiq¬ 
uitous in systems of coupled oscillators. Future works should 
attempt to extend this theory to the non pertubative regime of 
large line capacity and include dissipation effects which be¬ 
come important as voltage angle differences become large. 

We thank R. Delabays for useful discussions. This work 
was supported by the Swiss National Science Foundation. 

Appendix A: Jacobi matrices 
I. Properties of Jacobi matrices. 

In this section we list some of the properties of the eigen¬ 
values and eigenvectors of Jacobi matrices. For further details 
and proofs of these results see Ref. [38]. Consider the follow¬ 
ing positive semi-definite, symmetric, tridiagonal matrix with 


strictly 

negative subdiagonals 




-bi 

0 

0 ' 


-b\ 

02 

-b2 0 


J = 

0 



0 





~^n-l 


, 0 


0 -bn-\ 

Qn J 


Such matrices are also known as Jacobi matrices and have the 
property that their eigenvalues are distinct [38]. Furthermore, 
the principal minors of J (the k* principal minor J is the trun¬ 
cated version of J consisting of Jij with ij = 1,...,^ < n) 
satisfy the following recurrence relation 

Dk^iiA) = - A)Dk(A) - bjDt-iiA), (A2) 

where Di;(A) is the characteristic polynomial of the k* princi¬ 
pal minor of J (Do(A) - 1 and Di(A) - ai - A). In particu¬ 
lar Dn(A) is the characteristic polynomial of J which vanishes 
when A is equal to one of the eigenvalues di < A 2 < ■ ■ ■ < A„ 
of J. It can be shown [38] that the sequence 

(D„_i(d), D„_ 2 (d),..., £»i(d), £»o(4)] (A3) 

contains j - I sign changes when evaluated at the y* eigen¬ 
value A = Aj. 

If = (Mj''\ ..., is the eigenvector of J associated to 
Aj it is straightforward to show that the coefficients of sat¬ 
isfy a recurrence relation which is similar to that of Eq. (A2) 

« + Okul^^ - bkU^ll^ = AjU-l ^, (A4) 

bkU'll^ = (ok - Aj)u[^^ - bk-iu)ll^ , 
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for A: e {1,... ,n} and ^ , = 0. In fact one obtains that 

* ’ ’ ^ U n+i 

(A4) is fulfilled by 


u 


ij) 

k 


oc 


b\. ..bk-i ' 


(A5) 


Hence, given the sign property of the sequence {D(/lj)), it is 
clear that the components of the eigenvector will have j -1 
sign changes. 

We mention a last property [38] of the eigenvectors of J 
which is useful for our electrical model. Given two eigenvec¬ 
tors and we have 


+ aku[^'’ - bku[^l^ = , 

-bk-iuf_^ + akuf - bkull^ = Aiuf . 


(A6) 


Eliminating ak one obtains 


bkiU/^ “,fc+i ^k+i^k )^bk-i\u^. Mj,_j ) - (Ai Aj)u^ Uy, , 

(A7) 

which summed over k - 1,2,... ,l gives 




>) = (d,-d,.) j] 


U^k^U^k ■ 

k k 


(A8) 


A = 1 


2. Connection to the chain model. 

The matrix -M constructed for our chain model prior to the 
line addition is a Jacobi matrix with the additional property 
that since a, = 7>, -i- 7>,_i its first eigenvalue Ai is equal to zero 
and the corresponding eigenvector is Thus, 

applying Eq. (A8) to -M with j = I and i - 2 yields a relation 
between any two consecutive components of 



Additionally, according to the properties of Jacobi matrices, 

there will be only one sign change in the list of the coefficients 

of which therefore will have the form 

Lastly, the orthogonality relation between and implies 

that 


u(1)Tu(2) ^0 O ^ m(2) ^ Q ^ 

A=1 

Using Eq. (AlO) and the sign properties of the coefficients of 
suffices to see that Eq. (A9) leads to the conclusion that 
the components of are monotonously ordered (i.e. for any 
i > j we either have or 


Appendix B: Perturbation theory for a generic graph 


Given a generic network of N nodes and lines of capacity K/j, 
let (0,) and [0, ) respectively denote the solutions of the power 
flow equations (1) before and after the addition of a line of ca¬ 
pacity S <Si Ki j between nodes a and p. Assuming the 0,’s are 
small deviations of the unperturbed solution (Qi 0, -H 50, with 
|b0,| 1), we expand the power flow equations to leading 

order in 6 

0 = Kn cos(0/ - 0i){66i - 66i) i ^ a, p, 

0 = ^ Kia cos(0; - 6a)i66i - 66a) + 6 sin(0^ -6a) i - a, 

l~a, li=j3 

0 = ^ Kk/j COS(0/ - 6/})(66i - 66fj) + 6 sin(0„ - 0^) i =p. 

l~P, li^a 

(Bl) 

Eqs. (Bl) can be rewritten using the stability matrix, M [de¬ 
fined below Eq. (5)], of the system prior to the line addition 

M66 - 6 sin 0„^v, (B2) 

where 66 - {66\,..., 66^) and v is the N dimensional vector 
whose component is equal to v, = 6ka - 6ip. 

M, being a real symmetric matrix, is diagonalized by an 
orthogonal matrix T whose /* column is the & eigenvector, 
of M. Eurthermore, the U(l) symmetry of the power 
flow equations implies that one of the eigenvalues of M is null 
(di = 0 assocuated to u*'* = (1,..., 1)/ VA). Since M is sin¬ 
gular it cannot be inverted. However, Eq. (B2) can be solved 
for the (50,’s by using the Moore-Penrose pseudoinverse of M 
defined as 


M-' = T 






(B3) 


where T - ..., and the Ai’s are the eigenvalues of 

M. is such that = 1 - where 

u(i)u(i) is equal to the N x N matrix having 1/A for all it’s 
entries. Multiplying (B2) by yields 


(501 ' 

1 

~ A 

' 1 . 166 ,' 


f M~^ - M~^ \ 

66n , 

, 'Ll 66, , 

= 6 sin 6 aj 3 

M~^ - M~^ 


The difference of (50,’s between any two nodes is given by 

66^ - 66J = ekj = b sin dafi [{m;^ - M,:j) - - M^)] , 

(B5) 

where the term 2/^^/ drops due to the global rotational in¬ 
variance of the power flow solution. Einally, Eq. (B5) can be 
expressed in terms of the eigenvectors of M making use of 
(B3). This yields 

eij = 6 sin ^ - uf') AJ^, (B6) 

l>2 


In this section we extend the calculation of the leading or¬ 
der correction to the Lyapunov exponent A 2 resulting from the 
addition of a new line to the case of a generic electric network. 


for any node i connected to node j. Having established the 
correction to the power flow solution (B6), it is straight for¬ 
ward to compute the leading correction to the stability matrix. 











AM, and to obtain the correction of the Lyapunov exponent of a generic network is presented in Eq. (18). 
AA 2 - The final expression for AA 2 in the case 
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